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ABSTRACT 

We present an analysis of the molecular hydrogen absorption system at Zabs = 2.811 
in the spectrum of the blazar Q 0528—250. We demonstrate that the molecular cloud 
does not cover the background source completely. The partial coverage reveals itself 
as a residual flux in the bottom of saturated H 2 absorption lines. This amounts to 
about (2.22 ± 0.54)% of the continuum and does not depend on the wavelength. This 
value is small and it explains why this effect has not been detected in previous studies 
of this quasar spectrum. However, it is robustly detected and significantly higher than 
the zero flux level in the bottom of saturated lines of the Lya forest, (—0.21 ± 0.22) 
per cent. The presence of the residual flux could be caused by unresolved quasar 
multicomponents, by light scattered by dust, and/or by jet-cloud interaction. The H 2 
absorption system is very well described by a two-component model without inclusion 
of additional components when we take partial coverage into account. The derived total 
column densities in the H 2 absorption components A and B are log A(H 2 )[cm“^] = 
18.10±0.02 and 17.82 ± 0.02, respectively. HD molecules are present only in component 
B. Given the column density, log A(HD) = 13.33 ± 0.02, we find A(HD)/2A(H2) = 
(1.48 ±0.10) X 10 significantly lower than previous estimations. We argue that it is 
crucial to take into account the partial coverage effects in any analysis of H 2 bearing 
absorption systems, in particular when studying the physical state of high-redshift 
interstellar medium. 
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1 INTRODUCTION 

As a result of their cosmological distances quasars 
(QSOs) appear as point-like objects. Various studies have 
aimed to explore the detailed inner structure of quasars, 
which is unresolved even for low-redshift active galactic nu¬ 
clei (AGNs), because of their remote distances and sub¬ 
parsec scales of their emission regions. In the standard AGN 
paradigm the central region is divided into an accretion disk, 
a dusty-torus, a jet, a broad line region (BLR), and a narrow 
line region (NLR). Each of these regions contribute differ¬ 
ently to the AGN emission spectrum. 


Direct imaging of the spatial structure of AGNs is pos¬ 
sible with current instruments mainly probing longer scales. 


of the brightest AGNs (e.g. 

Jaffe et al.ll2004l: iTristram et al.l 

I 2 OI 4 I: iLopez-Gonzaea et al. 

I 2 OI 4 I) have revealed the exis- 


tence of a hot, parsec-scale disk that is surrounded by warm 
dust extended in the polar direction. In the optical band 
the geometry of the emission line region is investigated by 
indirect methods. Reverberation mapping establishes the re¬ 
lationship between the size and the luminosity of the BLR 
and yields a typical BLR size of Rblr ~ 0.2 pc ([Kaspi et al.l 
12007 : [^elouche V Danielll2012l ) for high redshift luminous 
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quasars. Differential microlensing allows f or a constraint on 
the accretion disk size < 3 x 1 0~^ pc (|Blackburne et aP 
1 20 111 : I Jimenez-Vicente et al.|[2012ll and for an est imation of 
the size of the BLR ~ 0.1 pc (jSluse et al.lboill l. The ob¬ 
servations of gamma-ray emi ssion constrain th e size of a jet 
constituent to a few parsecs (|Abdo et al.ll2010l L 

Another estimate of the size of the AGN emitting re¬ 
gions comes from constraints derived from covering factor 
analysis of intervening H 2 bearing clouds which happen to 
cover the background source only partially. Analysis of the 
partial coverage of Q 1 232+082 by a molecul ar hydrogen ab¬ 
sorption cloud allowed [Salashev et al.l (|201l| j to estimate the 
size of the CIV BLR, R^iv ~ PC. 

Molecular hydrogen absorption systems, a subset 
of damped Lya systems (DLAs) and sub-damped Lya 
systems (sub-DLAs), reveal diffuse and translucent in- 
terstellar clouds i n hig h redshift intervening galaxies 
(|Noterdaeme et aP I2OO8I L An analysis of H 2 absorption 
systems allows for examining the physical conditions of 
diffuse clouds in dista nt galaxies (jSrianand et al.l l2005l : 
iNoterdaeme et all l2007l j . ft has been shown that the gas 
is a part of the cold neutral medium with comparatively 
low kinetic temperature (T ~ 50 — 500 K) and high den¬ 
sities (riH > 10 cm“^), thus compact sizes (I < Ipc). 
Comparison of 21 cm and H 2 absorptions suggests that 
the H 2 absorption originates from a compact gas that 
probably contains only a small fraction of Hi measured 
along the line of sight (jSrianand et al.l l2012l j . These sys¬ 
tems are important instruments for the analysis of sev¬ 
eral cosmological problems, as follow s, (i) The discov¬ 
ery o f HD/H 2 clouds at high redshift (|Varshalovich et ^ 
l200lh provides an independent way to estimate the pri¬ 
mordial deuterium abundance (D/H) and therefore the 
relative baryon density of the Univ erse Qb which i s one 
of the key cosmolog ical parameters (jBalashev et al.l I2OIOI : 
llvanchik et al.l |2Q1 qI l (ii) The comparison of H 2 wave¬ 
lengths observed in QSO spe ctra with laboratory ones 
(just fo r this quasar Q 0528 — 250. | Varshalovich V LevshakovI 
(I 1993 IL ICowie &: Songaifl (I1995 Il IPotekhin et al.l (|l998l L 
lUbachs V Reinholdl (|2004l j. iKing et al.l (|201ll U allows us to 
test the possible cosmological variation of the proton-to- 
electron mass ratio = mpime- Because the de-composition 
of H 2 absorptions into several components is crucial for 
studies of the fundamental constant variability problem, 
we should pay attention to the partial coverage effect, ft 
is known that taking into account the partial coverage ef- 
fects, the physical mo del of the absorption system differs (see 
iBalashev et al.|[201lh . (iii) The interpretation of the relative 
populations of C i fine-structure excita t ion levels and CO 
rotat ional levels (jSrianand et al.l l200(]l : INoterdaeme et ahl 
l201l| j allows us to measure the temperature Tcmb {z) of the 
cosmic microwave background radiation at high redshift. 

Here, we argue that it is necessary to take into account 
the partial coverage of quasar emission regions by a com¬ 
pact intervening H 2 cloud in order to derive a robust fit 
of the absorption lines, ff this effect is not taken into ac¬ 
count, column densities can be underestimated by a factor 
of up to two orders of magnitude . The first case of s uch 
an analysis has been presented by IBalashev et al.l I 2 OIII for 
Q 1232+082 . The second case of par t ial co verage has been 
detected bv lAlbornoz Vasauez et aT ] (|2014| j for H 2 bearing 
cloud towards the quasar Q 0643—504. The third case of par¬ 


tial coverage for the H 2 cloud at Zabs = 2.811 in the spec¬ 
trum of Q 0528—2508 is presented in this study. We anal¬ 
yse a new spectrum and detect residual flux in the bottom 
of saturated H 2 lines (J = 0 and J = 1 levels). In case this 
flux is not taken into account, saturated lines yield large 
values and a multicomponent model is used instead (e.g. 
iKing et ^l201lh . 

The remainder of this paper is organized as follows. A 
brief description of the data is given in Section!^ The prin¬ 
ciples of partial coverage are described in Section[3l In Sec- 
tionUl we present the analysis of the H 2 absorption system, 
accounting for partial coverage . The HD molecular lines are 
explored in Section^ The results are discussed in Section[6l 
and we give a brief conclusion in Section[71 


2 DATA 

The molecular hydrogen was identified for the first 
time at high redshift in the very spectrum of Q 0528—2508 
(|Levshakov Varshalovichlll985l i . This quasar was observed 
many times, in particular during the period between 2001 
and 2009 using both spectroscopic arms of Ultraviolet and 
Visual Echelle Spectrograph (UVES) of the Very Large 
Telescope (VLT); f or a description of the instrument, see 
IPekker et al.l (|2000l i. The log of the observations used in our 
work is shown in Table [1] These observations relate to four 
programmes, three of which were carried out in 2001-2002: 
66.A-0594(A) (PI: Molaro), 68.A-0600(A) (PI: Ledoux), and 
68.A-0106(A) (PI: Petitjean). The instrument settings used 
during these observations were a 1-arcsec slit and 2x2 CCD 
pixel binning in both arms, resulting in a resolving power of 
R~45000 in the blue and R~43000 in the red. There was no 
ThAr lamp calibration attached to each of the exposures. An 
additional series of observations was performed in 2008-2009 
under programme 082.A-0087(A) (PI: Ubachs). The settings 
for that programme were a 0.8-arcsec slit in the blue arm and 
a 0.7-arcsec slit in the red. The 2x2 CCD pixel binning was 
also used at that time. This resulted in a resolving power of 
R^GOOOO in the blue and R~56000 in the red. Because the 
goal of that programme was to set a limit on the variation 
of /i, ThAr lamp calibrations were also taken immediately 
after each observation. 

The data presented in Table 1 were reduced using 
the UVES Common Pipeline Library (CPL) data reduction 
pipeline release 4.9.5 using the optimal extraction methocQ. 
The inter-order background (scattered light inside the in¬ 
strument) was carefully subtracted in both the flat-field 
frames and the science exposures. Linear spline interpolation 
was used to produce a two-dimensional background image, 
which was subsequently smoothed using an average boxcar. 
Eourth-order polynomials were used to find the dispersion 
solutions. However, the errors only reflect the calibration er¬ 
ror at the observed wavelengths of the ThAr lines used for 
wavelength calibration. All the spectra were corrected for 
the motion of the observatory around the barycentre of the 
Solar system and then converted to vacuum wavelengths. 
These spectra were interpolated into a common wavelength 


^ see the UVES pipeline user manual available for download at 

ftp: / / ftp. eso.org / pub / dfs / pipelines / uves / uves-pipeline-manual-22.8. pdf 
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Table 1. Log of the observations 


No. 

UT Date 

Program ID 

Exposure 

(sec) 

Slit 

(arcsec) 

1 

03.02.2001 

66.A-0594(A) 

1x5655 

1.0 

2 

04.02.2001 

66.A-0594(A) 

2x5655 

1.0 

3 

05.02.2001 

66.A-0594(A) 

1x5655 

1.0 

4 

07.02.2001 

66.A-0594(A) 

1x5655 

1.0 

5 

13.02.2001 

66.A-0594(A) 

1x5655 

1.0 

6 

13.03.2001 

66.A-0594(A) 

1x5655 

1.0 

7 

15.03.2001 

66.A-0594(A) 

1x5655 

1.0 

8 

17.10.2001 

68.A-0600(A) 

1x3600 

1.0 

9 

18.10.2001 

68.A-0600(A) 

2x3600 

1.0 

10 

08.01.2002 

68.A-0106(A) 

2x3600 

1.0 

11 

09.01.2002 

68.A-0106(A) 

2x3600 

1.0 

12 

10.01.2002 

68.A-0106(A) 

2x3600 

1.0 

13 

23.11.2008 

082.A-0087(A) 

2x2900 

0.8-0.7 

14 

25.11.2008 

082.A-0087(A) 

1x2900 

0.8-0.7 

15 

23.12.2008 

082.A-0087(A) 

4x2900 

0.8-0.7 

16 

25.01.2009 

082.A-0087(A) 

1x2900 

0.8-0.7 

17 

26.01.2009 

082.A-0087(A) 

1x2900 

0.8-0.7 

18 

26.02.2009 

082.A-0087(A) 

1x2900 

0.8-0.7 

Total 



108450 



array and generated the weighted-mean combined spectrum 
using the inverse squares of errors as weights. All the avail¬ 
able exposures were utilized to increase the signal-to-noise 
ratio up to ~60 per pixel in the wavelength range of the 
H 2 absorption lines (at z = 2.811). As shown below, this 
allows us to detect and study the effect of partial coverage 
in details. 


3 EFFECT OF PARTIAL COVERING 

Partial covering implies that only a part of the back¬ 
ground source is covered by the absorbing cloud. Mainly 
this can be a consequence of the absorbing cloud size being 
comparable to, or even smaller than, the projected extent of 
the background source. 

Partial covering is readily detectable in the spectrum of 
a background quasar if the cores of the saturated absorption 
lines do not reach the zero flux level. This indicates that a 
part of the radiation from the QSO passes by the cloud. The 
covering factor characterizing partial coverage is defined by 
the ratio. 


Tcloud 

Ttotal 


( 1 ) 


where iTioud is the flux that passes through the absorbing 
gas and Ttotai is the total flux. Therefore the measured flux 
in the spectrum, T(A), can be written as: 

F{X) = [Ttotai(A) — iTioud(A)] + Tcioud(A) exp[— t(A)], (2) 

here t(A) is the optical depth of the cloud at the wavelength 
A. The line flux residual (LFR) is the fraction of the QSO 
flux which is not covered by the cloud. These definitions are 



Figure 1. An illustration of the effect of partial coverage on the 
absorption-line profiles. Different panels show (a) a highly satu¬ 
rated line with total coverage, /c = 1, (b) a highly saturated line 
with partial coverage, fc = 0.8 and (c) a partially saturated line 
with the same partial coverage, fc = 0.8. The line flux residual 
(LFR) is the fraction of the QSO flux that is not intercepted by 
the cloud. It can be easily derived from the spectral analysis in 
case (b), but it is not straightforward to detect in case (c). 


illustrated in Fig.[T] The determination of the covering fac¬ 
tor is trivial in the case of highly saturated absorption lines 
(see Fig.[T]panels (a) and (b)) while for a partially saturated 
line (see Fig. [T] case (c)) the analysis requires a more sophis¬ 
ticated procedure. In this case it is necessary to use several 
absorption lines originating from the same levels but with 
different values of A/, which is the product of the oscillator 
strength, /, and the wavelength of t he transition, A . Such 
an analysis has been performed by Ivanchik et al.l (|2010l j 
and more precisely bv lBalas hev et ^ I (l201lh for the spec - 
trum of Q 1232+082, and bv lAlbornoz Vasauez et al.l (|2014| j 
for the spectrum of Q 0643 —504. A sim il ar sit uation was 
observed for HEOOO1—2340. I Jones et~aI1 (|2010l j have con¬ 
sidered the possibility of partial coverage of the BLR to ex¬ 
plain the observed Mgii equivalent widths. In contrast to the 
rare situations where partial covering occu rs from interven¬ 
ing systems (see also IPetitiean et al.lHoOOl j , partial covering 
is t ypical for absorption systems associated with quasars ( 
e.g. IPetit iean et Jll994 Rupke et al.ll2005l : iHamann et alJ 
l2ninl:lMiizahid et, al.lbniil. 

A failure to take into account the partial coverage effect 
in a spectroscopic analysis can lead to a significant underes¬ 
timation of the column density of an absorber. The system¬ 
atic bias (of column density) can exceed several orders of 
magnitude for saturated lines. As an example, consider an 
absorption line which consists of one component and has the 
high column density. The spectrum and the corresponding 
one-component model are shown in panel (a) of Fig. In 
panels (b), (c), (d) and (e) the same line is presented, but 
part of radiation from background source (10% for clarity) 
passes by a cloud. If we take into account the LFR, then the 
line can be properly fitted by a one component model and 
we can recover the high input column density (with given 
accuracy) and measure the LFR value (panel (e) of Fig. [2]). 
If the residual flux is not taken into account, then the one 
component model - Lorentzian (b) or Gaussian (c) profiles 
- is not adequate, the reason being that a one-component 
Voigt profile cannot describe the unsaturated bottom and 
far wings of the line simultaneously. Using additional com- 
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Figure 2. Panel (a) shows a single-component H 2 absorption line with log A/” = 18.2 and b = 4kms“^ and a fit by a one-component 
model (red line). The flux at the bottom of the line goes to zero. We add a residual flux of about 10 per cent (for clarity) to the data in 
panels (b), (c), (d) and (e). The flux at the bottom of the lines does not go to the zero level, and simultaneously the lines have Lorentzian 
wings. A simple one-component model cannot describe this line (cases (b) and (c)) and returns a large reduced x^- However, if we add 
a ne w unresolved com ponent to the model, the x^ will be significantly decreased (case (d)) but not the statistical criterion, the AlCC 
(see iKing et ^1201 ih , and the returned column density will be two order of magnitude too small. 


Table 2. Results of the previous analyses of the H 2 system at 
2 ; = 2.811 towards Q 0528—250. 


Year 

log Atot 

Ncomp. 

Resolution 

Ref. 

1985 

16.46T0.07 

1 


[1] 

1988 

18.0 

1 

10 000 

[2] 

1998 

16.77T0.09 

1 

10 000 

[3] 

2005 

2006 

18.22to\\r 

18.45T0.02 

2 

40 000 

[4] 

[5] 

2011 

16.56T0.02 

3 

45 000 

[6] 

2015 

18.28T0.02 

2 

45 000 

This work 


[1] ^evshako^^^^rshalovi^ (T9^); [2] ^oltz_e^^L (T98^; 
[3] ^rianan^^^etitie^ (T99^; [4] ^rianan^_e^^L (200^; 
[5] ^irkovic_e^^L (2006); [6] JCing_e^^L (201^ 


ponents, as shown in panel (d), a result with a satisfactory 
X^ is obtained. However, this solution is incorrect, because 
the resulting column density (logiV = 15.8) is much smaller 
than the input one (log A" = 18.2). To distinguish between 
cases (d) and (e) we propose a new method based on an 
analysis of several absorption lines with different oscillator 
strengths. The description of the method and application 
to the analysis of H 2 in Q 0528—250 are described in more 
detail later. 


4 MOLECULAR HYDROGEN 

Molecular hydrogen lines are detected in the spectrum of 
Q 0528—250 from the DLA at Zabs = 2.811. Column den¬ 
sity of neutra l hydrogen in this DLA s ystem is log A(Hi) = 
21.35 zb 0.07 (jNoterdaeme et al]l2008l ). H 2 lines correspond 
to transitions from rotational levels up to J = 5. To fit the 
molecular hydrogen lines, the spectrum has been normal¬ 
ized with a continuum constructed by fitting the selected 
continuum regions devoid of any absorptions with spline. 



Figure 3. Examples of H 2 line profiles corresponding to the tran¬ 
sitions from J = 0 — 5 levels. Two components are clearly seen in 
the J = 4 and 5 lines. We have found that two components are 
enough to obtain a satisfying x^- The origin of the velocity scale 
is taken at the redshift of the H 2 component, za = 2.81099. 


4.1 Number of components 

The profiles of the H 2 absorption lines have a com¬ 
plex structure that cannot be fitted with a single com¬ 
ponent. At least two components are clearly seen in the 
lines corresponding to J = 4, J = 5 rotational levels (see 
Fig [3}. Sinc e the first identification of the H 2 system in 
the quasar ([Levshakov fc Varshalovidilll985l i other studies 
have b een conducted, provid ing discordant results (see Ta¬ 
ble [^ . ISrianand et aU (|2005l ) used a two-component model, 
while liUng et al.l(|201lh pointed out that a three-component 
model is ve r y stro ngly preferred over two-component model. 
iKing et ^ (| 20 111 ) used a fitting procedure where they in¬ 
creased the number of components in the absorption sys¬ 
tem in order to minimize the corrected Akaike information 
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Figure 4. Some of the absorption lines of H 2 from the J = 0,1 levels detected in the spectrum of Q 0528—250 and the best Voigt profile 
fit (red line). Two components of the absorption system are shown by blue dashed lines. Only unblended saturated lines are present. 
The presence of prominent Lorentzian wings of the lines indicate high H 2 column density. The additional panels show a close-up of the 
bottom of the corresponding lines. The dashed black horizontal line represents the zero flux level. It is clearly seen that about 2 per cent 
of the LFR is present at the bottom of the saturated H 2 lines. The x-axes show the velocity offset from the centre of H 2 component 
at z = 2.81099. In the profile of H 2 line L9R1, some pixels in the bottom of the line are lower than 2 per cent level, which is probably 
because this line is blended with Lya forest absorption lines. 


criterion, AICC0 (|Sugiuralll978l : [King et al']l201lh . The re¬ 
sulting to tal column dens i ty dif fers by two order of magni¬ 
tude from ISrianand et al.l (|2005l ) . In that case a criterion for 
choosing preferred model is the consistent derived physical 
parameters of a cloud. It can be noted note tha t by using 
H 2 column density reported bv iKing et al.l (|201lh (based on 
the three-component model) a A/'(HD)/ 2 A/'(H 2 ) ratio is ob¬ 
tained, that is about an order of magnitude higher than the 
primordial one. Mea nwhile, the H 2 column density in two- 
component model of ISrianand et aP (|2005l ) gives a reason¬ 
able N(HD)/ 2 N(}l 2 ) ratio, which is c onsistent with the typ - 
ical values measured at high redshift (jBalashev et al.|[201Ql ) . 

A very important point for the choice of a reasonable 
absorption profile model in the case of Q 0528—250 is the 
presence of Lorentzian wings in J = 0, J = 1 line profiles 
(see Fig. mi. It is an indicator of the high H 2 column den¬ 
sity of absorption system with log N (H 2 ) > 18 which is con - 
sistent with the result reported by ISrianand et al.l (|2005l ) . 
However, in the new spectrum of Q 0528—250 the H 2 lines 
with prominent Lorentzian wings have some residual flux 


^ This statistical criteria allows for a choice of prefered model 
among several models with different numbers of fitting parame¬ 
ters. AICC = 7 where p is the number of fitting 

parameters, n is the number of spectral points, included in an 
analysis. 


at the bottom which significantly differs from the zero flux 
level (see Fig. and Fig. [5]). As a consequence, the fit by 
ISrianand et al. ~ (|2005l ) gives the large reduced x^ • It is prob¬ 
able that the signal-to-noise ratio in the previous spectrum 
was insufficient to detect the residual flux. 

The current study shows that the residual flux detected 
in the bottom of saturated H 2 lines in the spectrum of 
Q 0528—250 is the result of a partial coverage effect. In this 
case, the profiles of H 2 lines can be very well fitted by two- 
component model. Therefore, there is no need to increase 
th e number of com ponents in H 2 profiles (such as it is done 
bv iKing et al.ll201ll ) to explain the complex structure of lines 
(an example is given in Fig. HI). In the next three subsections, 
we provide evidence of the existence of the partial coverage 
in the spectrum of Q 0528—250. 


4.2 Zero-flux level correction 

To measure the LFR in H 2 absorption lines we need to derive 
the zero flux level in the spectrum. A non-zero flux in the 
core of a saturated absorption line can be the result of inac¬ 
curate determination in between spectral orders of scattered 
light inside the instrument. The zero flux level in the spec¬ 
trum can be estimated using the saturated Lya absorption 
lines which are numerous and almost uniformly distributed 
over the wavelength range where H 2 absorption lines are lo- 
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Figure 6. An illustration of the residual flux measurement pro¬ 
cedure. To estimate the residual flux we determine the median 
flux for pixels in a line profile that are in the range of Tier from 
the flux at the nearest pixel to the central wavelength Ac- This 
region is shown by purple shading. We refer to the width of this 
region as the width of the bottom of the line. 



cated. Lya lines are associated with intergalactic clouds that 
are larger than several kpc, thus it is most likely that they 
cover the background source completely. 

Wide flat bottom lines (AA > 1 A, i.e. Sv > 80 kms“^) 
were selected, which guarantees that lines are saturated and 
therefore measured fluxes in the bottom of the lines (LFRs) 
are the real zero flux level in the spectrum. These lines were 
selected in the spectral region 3500—4700 A. To estimate 
the residual flux at the bottom of the line the procedure 
illustrated in Fig. [6] was implemented. We selected several 
pixels in a line profile for which flux is within /c di Icr^, where 
fc is the flux at the centre of the line and ai is the error in 
pixel i. The residual flux was calculated as the median of 
the flux in the selected pixels. The width of the line bottom 
was estimated as the difference between the right most and 
left-most selected pixels from the line centre (see Fig. El). The 
top-left panel of Fig.[5]shows the LFR values obtained for the 
selected saturated Lya absorption lines. The average value is 
found to be (—0.21 zb 0.04)%. The standard deviation of the 
points is cr ~ 0.22%. The green stars show the LFR values 
measured at the bottom of the Lya and Ly/3 lines associated 
with the DLA-system at Zabs = 2.79 and with the Lya line 
of the second DLA-system at Zabs = 2.14. These lines are 
the most saturated lines for the spectrum. 


4.3 Partial coverage of H 2 absorption lines 

To estimate the residual flux in H 2 lines we selected J = 0,1 
lines without apparent blends. The LFR was measured by 
the same technique as for the Lya forest lines. The analy¬ 
sis found that the flux in the bottom of the lines is quite 
constant over large velocity range Sv ~ 10kms“^ (i.e. the 
dispersion of points is within the range of the average statis¬ 
tical error) that is wider than the full width half maximum 
(FWHM) of the UVES (6kms“^) and is comparable with 
the shift between centres of the two components. Therefore, 
the structure of H 2 system has no effect on the residual flux. 


The comparison of the obtained residual flux in the J = 0,1 
H 2 lines with the level of the residual flux in the Lya forest 
lines is shown in Fig|5] The filled squares represent the resid¬ 
ual flux in H 2 lines versus its location in the spectrum. We 
have estimated the line flux residual at the level 2.22zb0.54 % 
of the continuum, which significantly exceeds the zero-flux 
level. 

Non-zero residual flux at the bottom of saturated lines 
can also be the result of the convolution of the satu¬ 
rated lines with the instrumental function, or an imperfect 
data reduction, and/or of the blend of several unresolved 
unsaturated lines. First, however H 2 lines are wide with 
widths larger than the FWHM of the UVES spectrograph 
(6 km/s, i.e. ~ 0.08 A), so after convolution with the instru¬ 
ment function, the flux at the bottom of these lines must 
still go to zero. Secondly, the improper data reduction is 
not a viable explanation of the residual fluxes of saturated 
Lya lines are consistently equal to zero. Also we tested the 
dependence of residual flux on the width of the bottom of 
line for H 2 and for the Lya forest lines. This is shown in 
EigEU Several lines in the Lya forest were detected; these 
have similar widths as the H 2 lines and go to the zero level. 
The profiles of these lines are shown in Eig. lB2l 

Lastly, the residual flux can not be a result of the com¬ 
position of several unresolved unsaturated lines because the 
lines have Lorentzian wings. 


4 . 3.1 The Xf test 

An additional test which confirms the presence of the partial 
coverage effect is discussed here (and illustrated in Eig.[7|). 
Molecular hydrogen lines from J = 0,1 levels have very dif¬ 
ferent values of the product Xf (from 0.6 for LOPl to 36.0 
for WlQl). ft is known that the flux in the bottom of an 
absorption line decreases exponentially when Xf increases, 
F oc exp (—A/). The results of the calculation of this de¬ 
pendence for an absorption line in the spectrum with VLT 
resolution (EWHM = 6kms“^) are shown in Eig.[7|by the 
dashed lines. Eor simplicity, the modeled line profile consists 
of one component. The Doppler parameter b = 4kms“^ and 
the Damping width Tim is the same as for the H 2 line L2P1. 
Because the equivalent width of the line gets to the log¬ 
arithmic part of the curve of growth, we do not consider 
the different Doppler parameters. Two dashed curves (violet 
and grey) are calculated for the line with column densities 
log V =15 and 16. The curve is shifted from right to left as 
the column density N increases. Eor a higher column den¬ 
sity, the residual flux equals zero for a wide range of Xf 
(because the line is saturated) and differs from zero only for 
small values of A/ ^ 5 (the case of an unsaturated line). 
In a case of a multicomponent line (i.e. the blend of two or 
three lines with small column density log N = 15), the resid¬ 
ual flux in the bottom of the blended line would also behave 
exponentially (dashed-dot curves). However, the behavior of 
the residual fluxes in the bottom of the saturated H 2 lines in 
the spectrum of Q 0528—250 is quite different. The residual 
fluxes in the H 2 lines are shown by the filled blue and red 
squares. The squares do not follow the expected behavior, 
moreover the points scatter similarly around a median value 
in the whole range of Xf. 

Two models were applied to obtain a consistent and cor¬ 
rect fit to the H 2 lines. Model (i) considers the best fit of 
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Figure 5. Top left: Analysis of the zero flux level in the spectrum of Q 0528—250. The LFRs measured at the bottom of the saturated 
absorption lines in the Lyo; forest are plotted against the wavelength as fllled squares. The residual fluxes are expressed in per cent of 
the continuum. The green stars show the LFR measured in the Ly/3 and Lya absorption lines of the DLA system at = 2.79 and the 
Lya absorption from the DLA system at Zabs = 2.14. The average value of the zero flux level is shown by the red horizontal line. The 
scatter in the LFR at the bottom of the lines (i.e. noise in the spectrum) is found to be at the level of 0.22 per cent. Bottom left: a few 
examples of saturated Lyo: lines. The third bottom panel is the zoom of a part of the second panel. Top right: the residual fluxes at the 
bottom of the H2 absorption lines from rotational levels J = 0 (blue) and J = 1 (red) are plotted against wavelengths as fllled squares. 
The average value of the LFR in the H 2 lines is show n by the blue line. Blue c ircles and orange horizontal lines show the positions of 
QSO broad emission lines and their widths taken from IVanden Berk et al.l (|200ll L Bottom right: A few saturated H 2 absorption lines are 
present. The red dashed lines show the estimated residual flux in each H 2 line. The black line indicates the zero level. The first panel 
shows the zoom of a part of the second panel. It can be seen that the flux in the pixels located at the bottom of saturated H 2 lines 
systematically departs from zero. The shift is about 2 per cent of the continuum. 


the certain H 2 line (e.g. L2P1), which has a value of Xf = 
4.23 and LFR = 1.7% of the total flux. The line has a wide 
flat bottom ((^r^ ~ 10kms“^) and Lorentzian wings. To fit 
this line without the partial coverage, it is necessary to use 
several unsaturated components in the line profile. There¬ 
fore, this model describes the line profile well. Model (ii) 
takes into account the partial coverage, and the line L2P1 
can be fitted by the one-component model with high column 
density log N = 18.2 and LFR = 1.7% (the blue line in the 
right-hand top panel). This model also describes the line pro¬ 
file accurately. Using only one H 2 line we cannot determine 
the most probable model. However, if we consider several 
H 2 lines from the same J level which cover a wide range of 
A/, we will be able to discriminate it. It is known that the 
lines from the same J level correspond to the same physical 
region of a molecular cloud; therefore the lines are described 
by the same set of physical parameters (the column density 
N and Doppler parameter b). The difference between line 
profiles originating from the same J level is caused only by 
the different values of Xf for the lines. Therefore, the correct 
model of an absorption system must describe all measure¬ 
ments of residual flux in the bottom of lines from a given J 
level simultaneously. Using the best fit parameters for mod¬ 
els (i) and (ii) we have calculated the dependence of the 
residual flux on Xf (thick red and blue curves, respectively). 
The thick red line cannot describe all squares in the main 
panel simultaneously, whereas the blue thick line can. 


4.4 Voigt profile Fitting 

A Voigt profile fitting of the H 2 absorption lines was per¬ 
formed, taking into account the partial coverage. To describe 
a complex structure of line profiles we have divided the to¬ 
tal flux detected by an observer into two parts: from a main 
source and an additional one. In addition, we consider the 
H 2 system to be composed of two components (A and B) 
at redshifts za = 2.81099 and zb = 2.81112. The light from 
the main source (~ 98% of the total flux) is intercepted by 
the two H 2 components and does not produce any residual 
flux in the saturated lines. The light from the additional 
source that passes by H 2 clouds is not absorbed and there¬ 
fore produces a uniform residual flux in H 2 lines. The flux 
in the absorption line of two components A and B can be 
described as 

Fix) = (1 - fc)Ftotai + (3) 

where fc (in relative units) is the covering factor for H 2 lines. 

However, because the physical conditions in the A and 
B clouds (such us linear size, volume density, etc.) might 
be different, we can expect the covering factors of quasar 
emission regions by two H 2 clouds also to be different. In 
this case the construction of the H 2 line prohles is more 
complicated and we present an analysis of this case below 
(see Appendix lX)) . Here, it is important to note that taking 
into account two covering factors does not allow for a better 
ht to the H 2 lines (see the discussion in Appendix 

Then the absorption lines for each J level were de- 
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Figure 7. The residual flux at the bottom of the H 2 absorption lines is given versus the parameter A/ (where / is the oscillator strength 
of a line). Blue and red points correspond to the H 2 lines from the J = 0 and J = 1 levels. The theoretical flux at the bottom of the 
one-component line with a fixed column density N and Doppler parameter b is shown by the dashed curve. The violet and grey curves 
correspond to column densities logAr(H 2 ) = 15 and 16. It is seen that all points cannot be described by one curve simultaneously. The 
blue dash-dotted line represents the calculated residual flux at the bottom of a composite line, which consists of two components with 
logA/’(H 2 ) = 15 for both lines and a velocity separation of Sv = 2kms“^. The green dash-dotted line presents the same model for three 
components with the same column densities and velocity separations oi 6v = 2 and 4 kms“^, respectively (see text). An example is 
shown in the right-hand top panel. We fit a certain H 2 line (L2P1) by two models: one component with the LFR (blue line) and three 
components (red line). Both models give adequate fit of the line. But model without the LFR cannot describe flux in all lines together, 
as shown by the thick red line in the graph. In contrast, all points are described by a one-component model with residual flux (see the 
thick blue line). 


scribed using seven fitting parameters: za, zb, bA, be, Na, 
Nb, fc- We used uniform values of fc over the whole wave¬ 
length range, because the residual flux in the H 2 lines 
is almost independent of the wavelength (see Fig.[5|). The 
Doppler parameter b is a function of the rotational level 
J. To estimate fitting parameters Markov Chain Monte 
Carlo (MCMC) method was implemented, and to speed 
the convergence the Affin e invariant ensemble sampler by 
ICoodman fc Weax3 (|2010l ) was used. The main advantage 
of this searching algorithm is to better explore a parameter 
space and to avoid using the partial derivatives of the 
function which eases a number of numerical issues. This al¬ 
lowed for more reliable estimates of the fitting parameters 
in comparison with other algorithms. 

4.5 Fitting results 

The H 2 absorption system at z = 2.811 has more than 130 
absorption lines from J = 0 to J = 5 rotational levels. For 


the current analysis, we selected H 2 lines that are free of 
any obvious blends. The sample examined contains 99 lines. 
The best fit of H 2 lines are shown in Fig lCll -Fig lC5l ranked 
following wavelength positions. The best-fitting parameters 
are illustrated in Table O The reduced is 1.08 (the num¬ 
ber of fitting points ~ 1500). 

The H 2 absorption system is highly saturated. The total 
H 2 column densities are 18.10 zb 0.02 and 17.82 zb 0.02 for 
the A and B components consequently. This is consistent 
with the presence of the Lorentzian wings in the profiles 
of J = 0 and J = 1 levels. The obtained orto-to-para ratios 
are 2.7 zb 0.1 and 3.2 zb 0.3 for the A and B components, 
respectively. The corresponding kinetic temperatures of the 
H 2 clouds are Toi,a = 141 zb 6 K and Toi,b = 167 zb 13 K. 

Fig.[8]shows a comparison of the H 2 excitation diagrams 
for our measurements and those of previous work. The left- 
hand panel shows the result of the analysis performed by 
iKing et ] (1201 ih . the right-hand panel shows the result 
from the present work. Because the results presented here 
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Table 3. Column densities and Doppler parameters obtained 
from Voigt profile fitting of the H 2 system at Zabs = 2.811 to¬ 
wards Q 0528—250 after taking care of partial coverage. 


System 

J 

^'abs 

log A (cm ^) 

b (km s ^) 

A 

0 

2.8109950(20) 

17.50 


0.02 

2.66 

± 0.05 


1 

2.8109950(20) 

17.93 


0.01 

2.71 

± 0.05 


2 

2.8109952(5) 

16.87 


0.03 

2.75 

± 0.03 


3 

2.8109934(5) 

15.97 


0.07 

2.87 

± 0.07 


4 

2.8109938(8) 

14.18 


0.01 

4.79 

± 0.11 


5 

2.8109938(8) 

13.58 


0.02 

5.04 

± 0.39 

B 

0 

2.8111240(20) 

17.16 


0.03 

1.17 

± 0.06 


1 

2.8111230(20) 

17.67 


0.02 

1.14 

± 0.06 


2 

2.8111235(7) 

16.64 


0.03 

1.22 

± 0.02 


3 

2.8111238(6) 

16.24 


0.06 

1.25 

± 0.03 


4 

2.8111231(6) 

14.20 


0.01 

1.72 

± 0.09 


5 

2.8111231(6) 

13.60 


0.02 

2.38 

± 0.45 


are close to data reported bv ISrianand et aP (|2005l l. these 
are not shown here. We show excitation diagram s for two H 2 
compo nents at za = 2.81099 (2.811001 found bv iKing et al.l 
(l201lhl and zb = 2.8 1112. The third H 2 component found 
bv king et al.l (|201lh at zc = 2.8109346(11) is represented 
by the green stars. Although the ratio of J = 0 and J = 1 
levels of third component is the same as in other compo¬ 
nents, the excitation diagram is not physically realistic (the 
excitation temperatures T 02 and T 13 for the third compo¬ 
nent are negative). This might be the result of the incorrect 
model being used. It is seen, that the dis crepancy between 
our results and those of iKing et al ] ll 201 lh is larger only for 
the low J levels, where the influence of the partial coverage 
effect on the structure of the line profiles is significant. For 
high J levels, where column densities of H 2 are less, the re¬ 
sults agree in la. To sum up, the total column density of 
H 2 from the measurement is 18.284 ± 0.025, that is about 
two orders of ma gnitude larger than the value reported by 
iKing et, alJ feOllh 16.556 ± 0.024. 

It should be noted, that the same values of N and b 
parameters for all transitions of one J level were used. Us¬ 
ing the model with taking into account the partial coverage 
effect we obtain the reduced — 1 without an increase of 
the statistical errors of the spectrum. This is important, be- 
cause in the prev ious analysis of H 2 system in Q 0528—250 
iKing et~n ] (1201 ih noted that without artificially increasing 
the statistical er rors, the reduced was ^ 1 (see the cap¬ 
tion of tabled in iKing et ahllioill i. 

The value of the residual flux in H 2 lines is fitted as 
an independent parameter of an analysis. The best value 
is 2.40 ± 0.07% of the continuum, which agrees with the 
average value of the residual flux obtained from the analysis 
of J = 0,1 H 2 lines (see Section lT3)) . 


5 HD MOLECULES 

In the new spectrum (obtained by co-adding all previous and 
additional new observations) the molecular HD lines asso¬ 
ciated with the H 2 absorption system were detected (us- 


Table 4. Best fitting parameters for HD molecular lines in the 
spectrum of Q 0528—250. 


^abs 

log A (cm ^) 

b(kms ^) 

2.811121 ±0.000002 

13.33 ±0.02 

2.25 ± 0.53 

HD/ 2 H 2 

(1.48 ±0.10) X 10-5 



ing improved laboratory wavelengths, llvanov et al ]|^. 
The prese nce of HD lines in this system was first re¬ 
ported by iKing et ^ ] il2niih . The HD molecular lines are 
present only in the component B. Some of the HD lines are 
shown in Fig.[9l We have estimated the HD column density 
logiV(HD)B = 13.33 ± 0.02 by analysis of the two most 
prominent unblended absorption lines, L4-0R(0) and L 8 - 
0R(0). The results of Voigt profile fitting are presented in 
TablelH Other HD lines are highly blended and cannot be 
used in the analysis. The obtained column density is sig¬ 
nificantly less than that required to produce self-shielding, 
log A'(HD) ~ 15, thus we can set only a lower limit to the 
isotopic ratio D/H in the cloud. Because the total column 
density of H 2 in the component B is log A/'(H 2 ) = 17.85 ± 
0.02, we estimate D/H ^ A^(HD)/ 2 A^(H 2 ) = (1.48 ± 0.10) x 
10“^. The obtained HD column density is clos e to the re¬ 
sult, l og A'(HD) = 13.267 ± 0.072, reported bv I King et al.l 
(l20ljl . However, taking into account the significantly larger 
H 2 column density of the component B we have obtained a 
lower value of N(HD)/2N(H2) than the result bv iKing et al.l 
(I2011I1 . This limit is consistent with D/H ratio obtained 
from the analyse s of atomic species in quasar spectra (e.g. 
lOlive et al.ll 2012 l i. The comparison of this result with other 
A'(HD)/ 2 iV(H 2 ) measurements at high redshift is shown in 
Fig.Hni The N(HD)/2N(H2) in this system is consistent 
with other values and correspond to predictio ns of deuterium 
chemistry model s of diffuse ISM clouds (e.g. iBalashev et al.l 
l2016l:lLisztl[20lil . 

Note that in the component B we detect HD and Ci 
whereas in the component A these species are not present, 
despite the higher H 2 column density of this component in 
comparison with the component B. As for the component 
A we set an upper limit for HD and C i column densities 
which are logA'(HD)A ^ 13.1 and logA'(Ci)A ^ 12.0. The 
lack of HD and C i in the component A might be the result of 
higher local UV radiation (e.g. bright young stars near cloud 
A). It can destroy HD and ionize C i while the H 2 molecules 
self-shield against the local stellar radiation because to their 
high column density. 


6 DISCUSSION 

The presence of a residual flux at the bottom of the 
saturated lines of the H 2 system at Zabs = 2.811 towards 
Q 0528—250 can be interpreted according to the arguments 
described in the following subsections. 

6.1 Unresolved multicomponent quasars 

At redshift z ~ 2 , the internal structures in the emitting 
regions of of quasars with transverse dimensions < 3kpc 
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Figure 8. The excitation diagram for H 2 towards Q 0528—250. Here, Nj is the column density of the transition from the J level with 
gj degeneracy. Solid lines (blue, red and green) correspond to the excitation tem perature Tqi derived from J = 0 and J = 1 levels. The 
left-hand panel shows the result of the analysis performed by iKing et~ ] (l 2 onh . The excitation diagrams for different components of 
the H 2 system are shown by blue, red and green colors. The right-hand panel shows the results from our present work. Instead of three 
unsaturated components we use two components with higher H 2 column densities. 
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Figure 9. The line profiles of HD molecules in the absorption 
system at Za,bs = 2.811 in the spectrum of Q 0528—250. HD lines 
are detected in component B only. The fit to the WO-OR(O) line is 
inconsistent owing to a poor definition of the local continuum near 
the base of the Lya line of the second DLA system at z^bs = 2.14. 


are unresolved by UVES observations. For example, bi- 
nary quasars with separations of the order of ~ lOkpc 
(|Hennawi et al.|[2006l : IVivek et al.ll2009l l might remain unre¬ 
solved. Partial coverage can arise if Q 0528—250 has a com¬ 
plex multicomponent structure and if not all of the compo¬ 
nents are covered by the absorbing H 2 clouds. 

The availab le Very Long Baseline Array (VLBA) image 
of Q 0528—250 (|Kanekar et akllioPol : ISrianand et a'i]l20ld ) 
shows an unresolved component containing ^ 94 per cent of 
the t otal flux in the radio band (see table 6 of ISrianand et all 
l2012l b Another ~ 6 per cent is probably emitted by a diffuse 
component. However, it is intriguing to note the consistency 
of these numbers with our findings. The spatial size of ra- 



logiV(H^) 


Figure 10. The measurements of A/’(HD)/A/’(H 2 ) vs N(H 2 ) in 
absorption sys t ems a t high redshift. The data are taken from 
llvanchik et al ] (l2015h . The estimates of A/'( HD)/ A/'(H 2 ) toward 
Q 0528—250 from this work and iKing et al.l (1201 ih are shown by 
red and blue circles, accordingly. The blue horizontal strip repre- 
sents the ratio of atomic D i and H i measured in quasar spectra 
(|01ive et al.ll2012ll . The solid dashed line corresponds to constant 
column density of HD, log ^(HD) = 15. 


dip e mission core component is 65 x 380 pc (|Kanekar et al.l 
l2009l b that is significantly larger than the size of the H 2 
clouds. We have also looked at the images of PKS 0528—250 
in 13 and 4 cm taken as part of National Radio Astron¬ 
omy Observatories’ VLBA calibrators. The sources are un¬ 
resolved even at higher resolution achieved in these images. 
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6.2 Dust scattering 

Scattering by dust is characterized by a narrow radiation 
pattern. About 90 per cent of the total flux is scattered to- 
wards the ob server within less than 5-deg opening angle (e.g. 
IPrainel[20031 b The size of the scattering region can be much 
larger than the size of an H 2 cloud. This is why, even in 
the case of total coverage of the QSO by an H 2 cloud, the 
scattered radiation can be registered as a residual flux in the 
bottom of H 2 lines. The scattered flux by a dust-rich region 
of a DLA system can be estimated as 

= (4) 

ILqSO 

where F is the flux that passes through an H 2 cloud and is 
registered by the observer, tdla is the mean optical depth 
along the line of sight toward the quasar, Qqso and Qdla 
are the solid angles (measured by an observer on Earth) of 
the QSO emission region and the scattering region, accord¬ 
ingly. f^Qso is the solid angle of the QSO emission region 
measured by an observer at the position of the DLA and 
Qsc is the solid angle in which most of radiation is scattered 
by dust. Therefore, the residual flux in the bottom of H 2 
lines is LFR = Fsc/F. Assuming that a typical dimension 
of the scattering region is as large as 1 — lOkpc, LFR can 
be of the order of several percent. This is consistent with 
observations. The scattered radiation can produce a non¬ 
zero residual flux for metal lines also. However, this is not 
detected for Q 0528-250. 

6.3 Jet-cloud interaction 

The jet emission pointing toward us contributes to the con¬ 
tinuum in the optical band. The extended part of the jet 
interacts with the external ISM and can warm up a cloud 
which is distant from the central source. Similar situation 
has been detected in the deep radio-optical- X-ray observa¬ 
tions of the radio galaxy PKS B2152—699 by ([Worrall et al.l 
l2012l b who have reported the first high-resolution observa¬ 
tions of the radio jet in the direction of an optical emission¬ 
line high ionization cloud. The Hubble Space Telescope im¬ 
age shows not only emission in the r egion of the high-ionized 
clo ud - interpreted as ion ized gas bv lTadhunter et 51 (|l988l j 
and iFosburv et all (|l998l j - but also emission associated with 
a radio knot. The measured optic al flux density of thi s knot 
at 4.97 X 10^^ Hz is 2.5 ± 0.2 /iJy (|Worrall et al.ll20l1 b This 
is strikingly similar to the value of the residual flux in the V 
optical band which is 2.4 ± 0.1 /xJy (u sing the apparent mag¬ 
nitude rriY — 17.34 of Q 0528—250, IV5ron-Cettv VeronI 
(I 201011 1 . The variability of a jet emission (during from a few 
days up to a month) should be high and should also induce 
variability of the residual flux in the H 2 line. This could be 
tested by further observations. 


7 CONCLUSION 

We have performed a detailed analysis of the H 2 system 
at Zabs = 2.811 towards the quasar Q 0528—250. This is a 
well-known sy stem and the first detection of H 2 molecules at 
high redshift ([Levshakov Varshalovi^ll985l j . The system 
has been analysed in a variety of research, but the results 


obtained are not consistent (see Table [2|) . The discrepancies 
can be explained by the presence of a residual flux in the 
bottom of H 2 lines, which has not been considered in previ¬ 
ous research. We have derived the mean value of the residual 
flux as (2.22 ± 0.54) per cent of the continuum. This is sig¬ 
nificantly higher than the zero flux level, (—0.21 ±0.22) per 
cent, determined by analysis of the Lya forest lines. 

Taking into account the residual flux in the H 2 lines 
we have obtained a consistent fit of the H 2 system using 
a two-component model with high column densities. The 
derived total column densities of components A and B are 
log A'(H 2 ) = 18.10 ± 0.02 and 17.82 ± 0.02, respectively. 

We have performed the analysis of HD absorption lines 
detected only in component B. The estimated column den¬ 
sity is logA^(HD) = 13.33 ± 0.02 and thus we derive 
iV(HD)/2Ar(H2) = (1.48 ±0.13) x 10“^ This value is consis¬ 
tent with other measurements of iV(HD)/2A/’(H2) in quasar 
spectra at high redshift and can be considered as a lower 
limit of the primordial deuterium abundance (Balashev et 
al. 2010). 

Some interpretations for the presence of the residual 
flux are being offered here: (i) a multicomponent quasar; (ii) 
scattering by dust; (iii) a jet-cloud interaction. We favour the 
latter interpretation (iii). However, new optical and radio 
observations of Q 0528—250 are necessary to confirm this 
and reject the others. 

We argue that taking into account partial coverage ef¬ 
fects is crucial for any analysis of H 2 bearing absorption 
systems in particular when studying the physical state of 
high-redshift ISM. 
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APPENDIX A: TEST OF A MODEL WITH 
TWO ADDITIONAL SOURCES OF RESIDUAL 
FLUX 

In this section, we consider a model of absorption system, 
where the projected area over the illuminating source is dif¬ 
ferent for components A and B. In that case a structure of H 2 
line profiles is more complicated. To describe a flux detected 
by an observer, we divide a quasar emission region into three 
parts: the main source and two additional sources. The main 
source is covered by both H 2 clouds and does not produce 
the residual flux in absorption lines. Additional sources can 
produce the residual flux in two different ways (depending 
on the geometry of additional sources and H 2 clouds, see 
Fig.EU as follows. 

(i) One source (not covered by both clouds) produces 
the same residual flux in profiles of both H 2 systems. An¬ 
other source is covered by one cloud and produces residual 
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main source 



„/ 

additional sources 



absorption cloud A 


absorption cloud B 


Figure Al. An illustration of the possible explanation of the 
partial coverage effect for a case with independent sources of ra¬ 
diation (see text for detail). 


Table Al. A Comparison between models of the H 2 system with 
a different number of additional sources of radiation. A indicates 
the model with one additional source (which used in this work). B 
and C correspond to models (i) and (ii) described in see App. m 
p given the number of fitting parameters, AICC is the Akaike in¬ 
formation criterion (see Sect. liH . m and n indicate the intensity 
of additional sources relative to the total flux of a quasar. 


Model 

P 

AICC 

AAICC 

m, 10-2 

n, 10-2 

A 

31 

3096.0 

12.9 

2.4 ±0.1 

0 

B 

32 

3083.1 

0 

2.0 ±0.1 

4.0 ±0.3 

C 

32 

3331.1 

248.0 

2.4 ±0.1 

7.6 ±0.3 


flux only in profile of the second system: 

F{X)=mFtotai{X) + nFtotai{X)e-^^^^A , 

+ (1 - m - 

(ii) Each source produces the residual flux only for one 
system: 

F{X) =mFtotai{X)e-^^^^'> + nFtotai{X)e-^^^^'> + ^ ^ ^ 

(A2) 

+ (1 - m - n)Ftotal{X)e~'"^^^'*~'"^^^\ 

where m and n indicate intensities of additional sources in 
relative units. The results of Voigt profile fitting of H 2 lines 
using different models of the H 2 system are presented in 
TablelATI The values of AICC for models A and B are sig¬ 
nificantly lower than for model C, which could indicate the 
presence of one additional source of quasar radiation which 
uncovered by both H 2 clouds. Taking into account the sec¬ 
ond additional source (for system B) does not dramaticaly 
changed the AICC value and we have not found strong evi¬ 
dence for choosing a preferred model. Therefore, we use the 
simplest model (with one additional source) to analyse the 
H 2 system. 


APPENDIX B: THE EXPOSURE SHIFT 
ANALYSIS 

Because the final spectrum of Q 0528—250 is the co-addition 
of 27 exposures the non-zero residual flux at the bottom 
of saturated absorption lines (such as H 2 lines) could arise 
due to the average velocity shifts between exposures, up 
to 500 ms“^, and/ or intra-order velocity distortions, up 
to 15 00 ms~^ (e.g. I Whit more et al.l I 2 OIOI : iRahmani et al.! 
I2ni3li . Also, such a small LFR could b e the result of the dra - 
matic shift of even one exposure (e.g. iRahmani et al.ll20lj ). 
In that case, all narrow saturated absorption lines would 



Figure Bl. The dependence of residual fluxes in the bottom of 
lines on the widths of the bottom of lines. Blue points correspond 
to all saturated Lya forest lines with A > 3500 A and reached the 
zero-flux level within 5%. The yellow points mark Lya absorbers 
which have also Ly/3 absorption lines in the spectrum. The red 
points represent the residual fluxes in H 2 lines. Among the Lya 
forest lines there are lines with LFR < 1 % and the width of the 
bottom < 0.3 A. We marked them by green colour. Existence of 
such lines indicate that the residual fluxes in H 2 lines are not the 
instrumental effect (see text for detail). 


have the same non-zero residual fluxes, whereas the wide 
saturated absorption lines (like most of the Lya forest lines) 
will go to zero. Eig. IbT] shows the residual fluxes obtained 
in absorption lines that have the residual flux within 5 per 
cent of the continuum and located in the region 3500-4700 A 
in comparison to their widths at the bottom. The measure¬ 
ment procedure was described in Section lT^ The red filled 
squares correspond to H 2 lines, and the blue, green and yel¬ 
low filled squares correspond to the Lya forest lines. The 
points marked in green represent the sample of the narrow 
saturated Lya forest lines which look like the H 2 lines, but 
the residual flux at the bottom of the lines goes to zero. The 
profiles of these lines are shown in Fig. lB2l The existence 
of such lines in the spectrum of Q 0528—250 implies that 
the dramatic shift between exposures does not exist, other¬ 
wise the residual flux would be present in these lines too. 
There are also several Lya forest lines that have non-zero 
residual fluxes comparable or even larger than the average 
H 2 residual flux. By performing visual inspections of these 
lines, Ly/3 absorption lines for four of them (yellow squares 
in Fig. lBlD were found, due to the fact that corresponding 
Lya lines have redshifts large than 2.43. For systems with 
z lower than 2.43 the Ly^^ absorption lines are located in 
the Lyman-break region in the spectrum. The Fig |B3l shows 
the result of Voigt fitting analysis of these lines. All of them 
have small HI column densities (log N < 14.3) and large 
Doppler parameters (b > 25kms“^). 


APPENDIX C: THE BEST FIT OF H 2 LINES 

Here we present the figures of the best fit of the H 2 lines. 
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V, km/ s 

Figure B2. The profiles of the Lya forest lines corresponded to green filled squares in Fig. lBlI For comparison we show also two F12 
lines, LlO-OP(l) and L4-0R(0), in the two bottom right panels. The red horizontal line represents the zero flux level. 
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Figure B3. The profiles of the Lya forest lines corresponded to yellow filled squares in Fig. lBlI The top panels show Lya absorptions. 
The bottom panels show corresponded Ly/3 absorptions. The best Voigt profile fit is shown by red lines. 



Wavelengths, A 

Figure Cl. The best fit of the H 2 absorption lines from J = 0 level at Zabs = 2.811 toward Q 0528—250 (after incomporating the 
correction for partial coverage). Two components of absorption system are shown by blue dashed lines. The residuals between fit and 
data are shown in the top of each panel by dark blue line. Horizontal green lines represent ±3 a levels. 
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Wavelengths, A 


Figure C2. The best fit of the H 2 absorption lines from J = 1 level at Zabs = 2.811 toward Q 0528—250. Colors and lines are the same 
as in Fignre lCTl 
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3584,5 3585,0 3556,5 3554,0 3554,5 


Figure C3. The best fit of the H 2 absorption lines from J = 2 level at z^bs = 2.811 toward Q 0528—250. Colors and lines are the same 
as in Fignre lCTl 
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'A. The best fit of the H 2 absorption lines from J = 3 level at z^bs = 2.811 toward Q 0528—250. Colors and lines are the same 
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Wavelengths, A 


Figure C5. The best fit of the H 2 absorption lines from J = 5 level at Zabs = 2.811 toward Q 0528—250. Colors and lines are the same 
as in Fignre lCTl 
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Figure C6. The best fit of the H 2 absorption lines from J = 4 level at Zabs = 2.811 toward Q 0528—250. Colors and lines are the same 
as in Fignre lCTl 
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